Simple relativistic model of a finite-size particle 
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Abstract 

Soluble model of a relativistic particle describing a bag of matter with fixed 
radius held together in perfect balance by a self-consistent combination of 
three forces generated by electromagnetic and massive scalar and vector fields 
is presented. For realistic values of parameters the bag radius becomes that 
of a proton. 
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1. Introduction 



"I never satisfy myself until I can make a mechanical model of a thing" . These words of 
Lord Kelvin [|l| describe best my motivation to search for a simple relativistic mechanical 
model of an extended charged particle. Since over the years the concept of a mechanical 
model has undergone a substantial evolution, I do not hesitate to include also fields as 
building blocks for "mechanical" models. 

In this Letter I present a very simple model of a relativistic charged object (one may 
think of it as a proton or a nucleus) that is solvable in terms of elementary functions. The 
model consists of a swarm of particles endowed with three types of charges e, gs, and gy 
interacting in a self-consistent manner with three relativistic fields. The constituent particles 
are occupying a bounded region in space (a bag) and are described by a scalar phase-space 
distribution function /(r, p, t). The solution of the field equations in this model can also be 
given a hydrodynamic interpretation in which case the bag becomes a droplet of pressureless 
fluid interacting with the three fields. An unexpected result of this investigation is the 
appearance of a quantization condition from which the radius of the bag (or the droplet) is 
determined. 

Almost a century ago Poincare P argued that to counterbalance electrostatic repulsion 
inside a charged particle one must introduce cohesive forces — Poincare stresses. In order 
to comply with the requirements of relativity theory, these stresses must possess a dynamics 
of their own, they must in Pauli's words "depend on physical quantities which are causally 
determined by differential equations" [Q . Nowadays we know that a possible candidate for a 
"causally determined physical quantity" that will hold a charged particle together is a scalar 
field. For unlike the electrostatic forces generated by vector fields, the forces between like 
charges generated by scalar fields are attractive. To prevent collapse one needs also a short 
range repulsion and that in turn can be supplied by a massive vector field. That explains 
the choice of the main ingredients of my model. They are, of course, being used quite often 
in particle physics and in nuclear physics 0]. 

Several classical models of a charged particle have been invented in the past (see, for 
example, and also the reviews of the classical electron theory p|-p!l|), but in all of them 



the Poincare stresses were introduced ad hoc. The relativistic model of a classical charged 
particle with a finite, sharply defined radius presented here is the first, to my knowledge, in 
which the Poincare stresses incorporate some realistic elements. It can be used in two ways. 
First, with its help one may illustrate and clarify some old problems of relativistic theories 



of extended objects, discussed already by Abraham [|T^, Poincare [0], Lorentz |T^, and von 



Laue [T^JT^, that also more recently continued to cause controversy |[16|JT7[| . Second, with 



the proper choice of parameters, the model may serve as a zeroth order approximation in 
theories of nuclei and their high-energy collisions. For it is clear today that if we are to use 
in realistic applications a semi-classical model of a relativistic extended object it will not be 
to describe an electron as Poincare and others had tried in vain. We may still try, however, 
to describe in this manner nucleons and nuclei whose quantum features are less predominant 
since their Compton wave lengths are much smaller than their physical dimensions. 
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2. Description of the model 



The starting point of my construction is the following set of relativistic equations de- 
scribing the motion of a relativistic dust of particles interacting in a self-consistent manner 
with the three fields (c = 1) 

[(m - gs(j)){dt + V • V) + mF ■ ^p] /(r, p, t) = 0, (1) 

d,F>^^ = ef, (2) 

{n + ml)(f) = gsp, (3) 

d,G^'' + ml.W = gvf. (4) 

where V and dp denote the derivatives with respect to r and p, respectively. The laboratory- 
frame three-velocity v and the three-force F are related to the spatial parts of the four- 
velocity and the four-force f^, 

r = eF^'^u, - gs{d^ - u^u^^d^^ + gyG^^u,, (5) 

in the usual way, namely, f* = u^/u^ = p^/Ep, = /vP. The antisymmetric field tensors 
and G^y are constructed in the standard manner from the corresponding potential four- 
vectors and W^. The four-current and the scalar density p are defined in terms of the 
distribution function as follows 

3^{v,t)= |^pM/(r,p,t), (6) 

P(r,t)= J^mf{r,p,t). (7) 

The set of coupled equations (|lD-(§) is a generalization of Vlasov- Maxwell equations 
used in plasma physics. Despite the nonrelativistic appearance of Eq.(|I|), full relativistic 
covariance of the theory based on these equations can be proven in the same manner as we 
have done earlier for the pure Vlasov-Maxwell theory. 



Equations (|l|)-(^) may be viewed as classical if one assigns the parameters ms and my 
a dimension of inverse length and uses them only as measures of the range of forces due to 
scalar and vector fields. Planck's constant will make its appearance at the end in Eqs.(p9D 
and (^0|) when the parameters of the model are expressed in terms of the masses of mesons 
and their coupling constants. 



3. A static solution of the model 



A solution describing a bag of finite size is most easily obtained in the static case, when 
all constituent particles are at rest, 

/(r,p,t)=5(p)p(r), (8) 
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and the field equations reduce to the following simple set (in the static case there is 

just one scalar density p since Ep = m) 



- AAo = ep, (9) 
(-A + ml)<p = gsp, (10) 
{-A + ml)Wo = gvp. (11) 

In order to satisfy also the equation (|l]) for the distribution function (H), I assume the 
following equilibrium condition 

pF = -pV(eAo - gs(p + gyWo) = 0. (12) 

This condition simply means that the net force acting on each particle in the bag vanishes. 
Note that the equilibrium condition is not imposed everywhere in space — that would lead 
to a trivial solution — but only in those regions of space where matter is present. That is 
why the bag's boundary must be well defined. The equilibrium condition severely restricts 
possible solutions, but fortunately it does leave room for some interesting ones. Since the 
solutions of the equations (^)-(|T^ can always be scaled, I shall normalize them by imposing 
the following normalization condition 

'"dVp = l. (13) 

A simple, spherically symmetric solution of Eqs.(|9|)- ([TT|) is written below separately for 
the inside and the outside of the bag. 
Inside {r < R): 

P = /+-/-, (14) 
eAo = e'i^-^)~Vo, (15) 
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(16) 



9vWo = gv{ ,o , o - 1.2 , ^2 > (17) 



Outside {r > R): 



p = 0, (18) 

e^o - ^, (19) 

^750=^e~"^^('^-^), (20) 

9vWo = ^e--^^^~^\ (21) 
47rr 

where 65, 6y, and Vq are constants and f± are the following S-wave solutions of the Helmholtz 
equations. 
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= ^ sin(^_ 

4:71 r 

The wave vectors k±. are determined from a biquadratic equation obtained from the equihb- 
rium condition (|T^, 

4 ^ (23) 

where 

Q' = e'~gl + g'y, (24) 



D = B'^- ie'^Q'^mlm^, (25) 

and 

5 = (^/l - e2)m2^ - {gl. + e^)ml. (26) 

The remaining parameters bs,bv, d±, the depth of the potential well Vq, and the radius of the 
bag R are determined from six continuity conditions at the bag's boundary. A combination 
of these conditions gives the following quantization condition for R, very similar to those 
arising in wave mechanics, 

T{ms)=T{mv), (27) 

where 

T(m ) = ^l + '^x k- + mx tanjk.R) 

kl +ni\k^ + mx tan(A;+i?) ' ^ ^ 

Eq.(p7D has infinitely many solutions for R, but only the lowest one is physically accept- 
able because all higher ones do not lead to a positive density p. Once the radius of the 
bag is determined, the remaining parameters can be calculated from explicit formulas. In 
particular, 

- 4^ TU-klt^/kl' 

where t± = tan{k±R) — k±R and T is the common value of T{ms) and T{mv). 

The solutions of field equations obtained in this way have several general features that 
are worth noting. The bag has a sharp edge — the value of the density at the surface of the 
bag is always finite, 

^^^^ = 

It is seen from this formula that for a solution to exist the effective repulsion must be stronger 
than attraction {Q^ > 0). 
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4. Hydrodynamic description 



My bag model may easily be converted into a droplet model by replacing the description 
in terms of the distribution function by a hydrodynamic description in terms of a density 
scalar field p and a four- velocity field that characterize the state of the fluid. The starting 
point of the hydrodynamic description is the following set of relativistic equations 



d^ipu^) = 0, (31) 

(m - gs<p)pu''d,u>' = epF'^'^u, - gsp{dy - u^u-'d^cP) + gvpG^^'u,, (32) 

d^F^-" = epu-", (33) 

(□ + m|)0 = gsp, (34) 

d^G"^ + mlW = gvpu\ (35) 



These equations describe matter modelled by a pressureless fluid of density p moving with 
four- velocity u^. The fluid is endowed with the three types of charges e, gs, and gv and is 
interacting in a self-consistent manner with the three relativistic fields. For a fluid without 
pressure the static solution remains the same as in the bag model. One may also add the 
pressure term and an equation of state and seek numerical solutions of the same general 
nature — with a sharply defined radius. 

The system described by the equations (|3lD -(^5D is conservative (and that is also true for 
the system described by generalized Vlasov- Maxwell equations (|l])-(^)) since the equations 
of motion guarantee that the total energy-momentum tensor T^'^ of the system, 

T'^- = (m - gs(p)pu^u'' + F^^F^' + ^g'^'F^pF^^ + dyd'^^P - g^'''\{dx<pd^<p - m|02) 

+ G^^G/ + mlW^'W'' + g^^'ilGxpG^P - \mlWxW^), (36) 

is conserved, d^T^^" = 0. 

The droplet (or the bag) described here is energetically stable. The combined energy of 
the three fields, owing to the equilibrium condition, reduces to half of the interaction energy 
with the scalar field (a version of the virial theorem) and that leads to the following value 
of the total energy of the system obtained from (^) by integrating the energy density 

Etot = J d?r{m - gs<p)p + \ j dh[{V A^f + (V^)^ + m|02 + {yw^f + mlW^] 
= J (fr{m - gs<p)p + \ j rfV[-AoAAo - + m|02 - Wo^^^ + mlW^] 
= m + \j d^r[eAQp - gs(j)p + gyWop] = m- Vo/2. (37) 

This formula conveys an important message concerning equilibrium configurations in local 
field theories. In order to calculate the total energy of the matter-field configuration in the 
whole space it is enough to know the field inside the bag. Since we know the relativistic 
transformation properties of the fields describing the model, a static solution may easily be 
boosted to give a model of a moving particle. The infamous factor of 4/3 discovered by 
J.J.Thompson [T^ that was plaguing all naive electron models does not appear here. It can 
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be checked by an exphcit calculation or inferred from a generalized virial theorem (cf., for 
example, 



that my solution satisfies the von Laue's condition [|15 

d^T'Hr) =0, 



(38) 



that guarantees proper transformation properties of the energy-momentum vector of the 
particle. 



5. A tentative connection with reality 



Now I come to the problem of assigning definite values to the parameters of the model 
in order to see if the bag can be made to resemble a nucleon. I have not done any elaborate 
parameter fitting, since that will definitely require the introduction of several scalar and 
vector mesons as is done in modern mean-field theories of nuclear structure |Q. I have just 
taken as my input the same values of the masses of the cx-meson and the cu-meson and the 
values of the coupling constants that are used in the simplest version (QHD-I) of the mean 
field theory of nuclear matter (Ref. [§l,p. 125), 

ms = 550 MeV /he, mv= 783 MeV /he, (39) 
gl = 91.64 he, gl= 136.2 he. (40) 

The quantization condition gives then the following values for the radius of the bag and 
the depth of the binding potential 

R = 1.05 fm, Vo = 31.42 MeV. (41) 

The values of the proton mean radius and the binding energy calculated for these values are 

rp = 0.714 fm, Eb = 15.71 MeV, (42) 

and they are to be compared with the experimentally measured radius 0.862 fm and the bulk 
binding energy per nucleon in nuclear matter 15.75 MeV. I am very far from advocating the 
use of this model in its present rudimentary form as a realistic model of nuclear structure 
and I am presenting these numbers only to show that they are not completely out of touch 
with reality. 

To illustrate my results I show in Fig. |l] separately the curves for the three potentials, 
in Fig. ^ the combined effective potential well, and in Fig. |^ the distribution of matter in 
the bag. All these graphs correspond to the choice of parameters given by the formulas ( |39D 
and ®. 

The potential well shown Fig. |^ may be used as an initial step in constructing a relativistic 
version of the nuclear shell model in which the Dirac equation would be used to describe 
matter. More ambitiously, one may use the understanding reached with my simple model 
(in particular, the quantization conditions for the radius) to solve more elaborate models in 
which matter is also described at the field theoretic level. 



7 



Acknowledgment 

I would like to acknowledge the hospitality extended to me by Walter Greiner at the 
University of Frankfurt. 



8 



REFERENCES 



[1] Lord Kelvin, Notes of Lectures on Molecular Dynamics and the Wave Theory of Light, 
(Johns Hopkins University, Baltimore, 1884), p. 270. 

[2] H.Poincare, Rendiconti del Circolo Matematica di Palermo 21 (1906) 129 (English trans- 
lation H.M.Schwartz Am.J.Phys. 39 (1971) 1287; 40 (1972) 862; 40 (1972) 1282). 

[3] W.Pauli, Theory of Relativity, (Pergamon, London, 1958), p. 186. 

[4] B.D.Serot and J.D.Walecka, The Relativistic Nuclear Many-Body Problem, Advances in 
Nuclear Physics Vol.16, edited by J.W.Negele and E.Vogt, (Plenum, New York, 1986). 
[5] P.A.M.Dirac, Proc.Roy.Soc. (London) 268 (1962) 57. 
[6] J.Schwinger, Found.Phys. 13 (1983) 373. 
[7] I.Bialynicki-Birula, Phys.Rev.D 28 (1983) 2114. 

[8] F.Rohrlich, Classical Charged Particles, (Addison- Wesley, Reading, 1965). 
[9] A.LMiller, Albert Einstein's Special Theory of Relativity, (Addison- Wesley, Reading, 
1981). 

[10] P.Pearle, Classical Electron Models in Electromagnetism, Paths to Research, edited by 

D.Teplitz, (Plenum, New York, 1982). 
[11] A.D.Yaghjian, Relativistic Dynamics of a Charged Sphere, Lecture Notes in Physics 

Mil, (Springer, Berlin, 1992). 
[12] M.Abraham, Phys.Zeit. 5 (1904) 576. 

[13] H.A.Lorentz, The Theory of Electrons, (Teubner, Leipzig, 1909), p. 213. 

[14] M.von Laue, Ann.Phys. (Leipzig) 35 (1911) 524. 

[15] M.von Laue, Das Relativitatsprinzip, (Vieweg, Braunschweig, 1911). 

[16] T.H.Boyer, Phys.Rev.D 25 (1982) 3246. 

[17] F.Rohrlich, Phys.Rev.D 25 (1982) 3251. 

[18] I.Bialynicki-Birula, J.C.Hubbard, and L.A.Turski, Physica 128A (1984) 504. 
[19] J.J.Thomson, Phil.Mag. 11 (1881) 229. 



9 



FIGURES 



FIG. 1. Separate graphs of the three self-consistent potentials plotted for the values of the 
parameters given by Eq.(|39|) and Eq.(|40|). Two repulsive potentials, long range electrostatic (owing 
to the smallness of the fine structure constant barely different from zero on this scale) and short 
range mesonic, are due to the vector fields. The attractive potential is due to the scalar field. 

FIG. 2. Total effective potential due to the combined effect of the three fields producing a 
potential well that keeps the constituents inside the bag. 

FIG. 3. The density of matter p for a spherically symmetric bag plotted as a function of r. 
Note a discontinuity of the density at the bag's boundary. 
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